Limitations of lymphoblastoid cell lines for establishing genetic reference datasets in the immunoglobulin loci.

Lymphoblastoid cell lines (LCLs) have been critical to establishing genetic resources for biomedical science. They have been used extensively to study human genetic diversity, genome function, and inform the development of tools and methodologies for augmenting disease genetics research. While the validity of variant callsets from LCLs has been demonstrated for most of the genome, previous work has shown that DNA extracted from LCLs is modified by V(D)J recombination within the immunoglobulin (IG) loci, regions that harbor antibody genes critical to immune system function. However, the impacts of V(D)J on short read sequencing data generated from LCLs has not been extensively investigated. In this study, we used LCL-derived short read sequencing data from the 1000 Genomes Project (n = 2,504) to identify signatures of V(D)J recombination. Our analyses revealed sample-level impacts of V(D)J recombination that varied depending on the degree of inferred monoclonality. We showed that V(D)J associated somatic deletions impacted genotyping accuracy, leading to adulterated population-level estimates of allele frequency and linkage disequilibrium. These findings illuminate limitations of using LCLs and short read data for building genetic resources in the IG loci, with implications for interpreting previous disease association studies in these regions.

Cutting Edge: Neutralizing Public Antibody Responses Are an Ancient Form of Defense Conserved in Fish and Mammals.

The repertoire of Abs is generated by genomic rearrangements during B cell differentiation. Although V(D)J rearrangements lead to repertoires mostly different between individuals, recent studies have shown that they contain a substantial fraction of overrepresented and shared "public" clones. We previously reported a strong public IgHµ clonotypic response against the rhabdovirus viral hemorrhagic septicemia virus in a teleost fish. In this study, we identified an IgL chain associated with this public response that allowed us to characterize its functionality. We show that this public Ab response has a potent neutralizing capacity that is typically associated with host protection during rhabdovirus infections. We also demonstrate that the public response is not restricted to a particular trout isogenic line but expressed in multiple genetic backgrounds and may be used as a marker of successful vaccination. Our work reveals that public B cell responses producing generic Abs constitute a mechanism of protection against infection conserved across vertebrates.

Large-scale analysis of 2,152 Ig-seq datasets reveals key features of B cell biology and the antibody repertoire.

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists.

An in vivo method for diversifying the functions of therapeutic antibodies.

V(D)J recombination generates mature B cells that express huge repertoires of primary antibodies as diverse immunoglobulin (Ig) heavy chain (IgH) and light chain (IgL) of their B cell antigen receptors (BCRs). Cognate antigen binding to BCR variable region domains activates B cells into the germinal center (GC) reaction in which somatic hypermutation (SHM) modifies primary variable region-encoding sequences, with subsequent selection for mutations that improve antigen-binding affinity, ultimately leading to antibody affinity maturation. Based on these principles, we developed a humanized mouse model approach to diversify an anti-PD1 therapeutic antibody and allow isolation of variants with novel properties. In this approach, component Ig gene segments of the anti-PD1 antibody underwent de novo V(D)J recombination to diversify the anti-PD1 antibody in the primary antibody repertoire in the mouse models. Immunization of these mouse models further modified the anti-PD1 antibodies through SHM. Known anti-PD1 antibodies block interaction of PD1 with its ligands to alleviate PD1-mediated T cell suppression, thereby boosting antitumor T cell responses. By diversifying one such anti-PD1 antibody, we derived many anti-PD1 antibodies, including anti-PD1 antibodies with the opposite activity of enhancing PD1/ligand interaction. Such antibodies theoretically might suppress deleterious T cell activities in autoimmune diseases. The approach we describe should be generally applicable for diversifying other therapeutic antibodies.

Single-cell analysis of human B cell maturation predicts how antibody class switching shapes selection dynamics.

Protective humoral memory forms in secondary lymphoid organs where B cells undergo affinity maturation and differentiation into memory or plasma cells. Here, we provide a comprehensive roadmap of human B cell maturation with single-cell transcriptomics matched with bulk and single-cell antibody repertoires to define gene expression, antibody repertoires, and clonal sharing of B cell states at single-cell resolution, including memory B cell heterogeneity that reflects diverse functional and signaling states. We reconstruct gene expression dynamics during B cell activation to reveal a pre-germinal center state primed to undergo class switch recombination and dissect how antibody class-dependent gene expression in germinal center and memory B cells is linked with a distinct transcriptional wiring with potential to influence their fate and function. Our analyses reveal the dynamic cellular states that shape human B cell-mediated immunity and highlight how antibody isotype may play a role during their antibody-based selection.

FATC Domain Deletion Compromises ATM Protein Stability, Blocks Lymphocyte Development, and Promotes Lymphomagenesis.

Ataxia-telangiectasia mutated (ATM) kinase is a master regulator of the DNA damage response, and loss of ATM leads to primary immunodeficiency and greatly increased risk for lymphoid malignancies. The FATC domain is conserved in phosphatidylinositol-3-kinase-related protein kinases (PIKKs). Truncation mutation in the FATC domain (R3047X) selectively compromised reactive oxygen species-induced ATM activation in cell-free assays. In this article, we show that in mouse models, knock-in ATM-R3057X mutation (Atm⁠ RX ⁠, corresponding to R3047X in human ATM) severely compromises ATM protein stability and causes T cell developmental defects, B cell Ig class-switch recombination defects, and infertility resembling ATM-null. The residual ATM-R3057X protein retains minimal yet functional measurable DNA damage-induced checkpoint activation and significantly delays lymphomagenesis in Atm⁠ RX/RX ⁠ mice compared with Atm⁠ -/- ⁠. Together, these results support a physiological role of the FATC domain in ATM protein stability and show that the presence of minimal residual ATM-R3057X protein can prevent growth retardation and delay tumorigenesis without restoring lymphocyte development and fertility.

Droplet-based mRNA sequencing of fixed and permeabilized cells by CLInt-seq allows for antigen-specific TCR cloning.

T cell receptors (TCRs) are generated by somatic recombination of V/D/J segments to produce up to 1015 unique sequences. Highly sensitive and specific techniques are required to isolate and identify the rare TCR sequences that respond to antigens of interest. Here, we describe the use of mRNA sequencing via cross-linker regulated intracellular phenotype (CLInt-Seq) for efficient recovery of antigen-specific TCRs in cells stained for combinations of intracellular proteins such as cytokines or transcription factors. This method enables high-throughput identification and isolation of low-frequency TCRs specific for any antigen. As a proof of principle, intracellular staining for TNFα and IFNγ identified cytomegalovirus (CMV)- and Epstein-Barr virus (EBV)-reactive TCRs with efficiencies similar to state-of-the-art peptide-MHC multimer methodology. In a separate experiment, regulatory T cells were profiled based on intracellular FOXP3 staining, demonstrating the ability to examine phenotypes based on transcription factors. We further optimized the intracellular staining conditions to use a chemically cleavable primary amine cross-linker compatible with current single-cell sequencing technology. CLInt-Seq for TNFα and IFNγ performed similarly to isolation with multimer staining for EBV-reactive TCRs. We anticipate CLInt-Seq will enable droplet-based single-cell mRNA analysis from any tissue where minor populations need to be isolated by intracellular markers.

Genomic analysis of a second rainbow trout line (Arlee) leads to an extended description of the IGH VDJ gene repertoire.

High-throughput sequencing technologies brought a renewed interest for immune repertoires. Fish Ab and B cell repertoires are no exception, and their comprehensive analysis can both provide new insights into poorly understood immune mechanisms, and identify markers of protection after vaccination. However, the lack of genomic description and standardized nomenclature of IG genes hampers accurate annotation of Ig mRNA deep sequencing data. Complete genome sequences of Atlantic salmon and rainbow trout (Swanson line) recently allowed us to establish a comprehensive and coherent annotation of Salmonid IGH genes following IMGT standards. Here we analyzed the IGHV, D, and J genes from the newly released genome of a second rainbow trout line (Arlee). We confirmed the validity of salmonid IGHV subgroups, and extended the description of the rainbow trout IGH gene repertoire with novel sequences, while keeping nomenclature continuity. This work provides an important resource for annotation of high-throughput Ab repertoire sequencing data.

Regulation of T-cell Receptor Gene Expression by Three-Dimensional Locus Conformation and Enhancer Function.

The adaptive immune response in vertebrates depends on the expression of antigen-specific receptors in lymphocytes. T-cell receptor (TCR) gene expression is exquisitely regulated during thymocyte development to drive the generation of αß and γδ T lymphocytes. The TCRα, TCRß, TCRγ, and TCRδ genes exist in two different configurations, unrearranged and rearranged. A correctly rearranged configuration is required for expression of a functional TCR chain. TCRs can take the form of one of three possible heterodimers, pre-TCR, TCRαß, or TCRγδ which drive thymocyte maturation into αß or γδ T lymphocytes. To pass from an unrearranged to a rearranged configuration, global and local three dimensional (3D) chromatin changes must occur during thymocyte development to regulate gene segment accessibility for V(D)J recombination. During this process, enhancers play a critical role by modifying the chromatin conformation and triggering noncoding germline transcription that promotes the recruitment of the recombination machinery. The different signaling that thymocytes receive during their development controls enhancer activity. Here, we summarize the dynamics of long-distance interactions established through chromatin regulatory elements that drive transcription and V(D)J recombination and how different signaling pathways are orchestrated to regulate the activity of enhancers to precisely control TCR gene expression during T-cell maturation.

T cell receptor repertoires of ex vivo-expanded tumor-infiltrating lymphocytes from breast cancer patients.

A higher level of tumor-infiltrating lymphocytes (TILs) is associated with better prognosis in breast cancer patients. Adoptive transfer of lymphocytes coupled with conventional therapies has appealed to many clinicians and investigators as an effective treatment strategy for cancer patients, which necessitates efficient activation and expansion of cytotoxic T lymphocytes precisely targeting cancer cells. To comprehensively understand composition of TILs and to provide a grounding in adoptive T cell therapy, we analyzed the T cell receptor (TCR) repertoires in ex vivo-expanded TILs from nine breast cancer patients via next-generation sequencing. For the three of them, TCR repertoires of TILs gathered after the initial culture during 2 weeks were additionally analyzed and compared to those of TILs that underwent ex vivo rapid expansion procedure (REP). Diversity of TCR repertoire was variable among the patients. V/J segment usage in the clonotypes was similar among patients, with variable distribution of read counts for each V/J segment. The top 50% of most frequently observed VJ combinations was present in > 80% of the total clonotypes. Compared with TCGA data, the samples contained a similar amount of recurrent CDR3 sequences, but clonotype expansion was variable among the samples. In terms of clinicopathologic factor, presence of in vitro reactivity among triple-negative breast cancer cases seemed to be related to lower Shannon's index, but p value was not statistically significant. In addition, the proportion of CD45RO+ cells out of CD8+ T cells were negatively correlated with Shannon's diversity index for both TCRα and TCRß chains (p = 0.010) via Spearman test. In this study, we identified a heterogeneous pattern of expanded T cell clones and stable usage of V/J segments in ex vivo-expanded TILs from breast cancer patients. Further large-scale studies are requisite to elucidate the clinical significance of TCR repertoires.

Templated insertions at VD and DJ junctions create unique B-cell receptors in the healthy B-cell repertoire.

RAG complexes recognise (cryptic) RSS sites both in and outside immunoglobulin sites. Excision circles may be reinserted into V(D)J rearrangements as long templated insertions to diversify the adaptive immune repertoire. We show that such VDJ with templated insertions are incidentally found in the repertoire of healthy donors.

Genetic interaction between the non-homologous end-joining factors during B and T lymphocyte development: In vivo mouse models.

Non-homologous end joining (NHEJ) is the main DNA repair mechanism for the repair of double-strand breaks (DSBs) throughout the course of the cell cycle. DSBs are generated in developing B and T lymphocytes during V(D)J recombination to increase the repertoire of B and T cell receptors. DSBs are also generated during the class switch recombination (CSR) process in mature B lymphocytes, providing distinct effector functions of antibody heavy chain constant regions. Thus, NHEJ is important for both V(D)J recombination and CSR. NHEJ comprises core Ku70 and Ku80 subunits that form the Ku heterodimer, which binds DSBs and promotes the recruitment of accessory factors (e.g., DNA-PKcs, Artemis, PAXX, MRI) and downstream core factors (XLF, Lig4 and XRCC4). In recent decades, new NHEJ proteins have been reported, increasing complexity of this molecular pathway. Numerous in vivo mouse models have been generated and characterized to identify the interplay of NHEJ factors and their role in development of adaptive immune system. This review summarizes the currently available mouse models lacking one or several NHEJ factors, with a particular focus on early B cell development. We also underline genetic interactions and redundancy in the NHEJ pathway in mice.

Monitoring minimal residual disease in the bone marrow using next generation sequencing.

Achieving minimal residual disease (MRD) negativity in the bone marrow is one of the strongest prognostic factors in multiple myeloma. Consequently, MRD testing is routinely performed in clinical trials and moving towards standard of care. This review focuses on the role of next generation sequencing (NGS) of tumor-specific immunoglobulin V(D)J sequences for MRD tracking. The immunoglobulin variable regions are ideal targets for tracking, because every tumor cell shares an identical gene sequence, which is stable over time and generally distinct from the immunoglobulin sequences of normal B-cells. Several excellent assays for NGS-based MRD testing are available, both commercial and community-based, including one that is FDA-approved. These assays can achieve the gold standard analytical sensitivity of one tumor cell per million (10-6), requiring a minimum input of 3 million bone marrow cells. On-going clinical trials will outline how MRD testing should be used to inform dynamic risk-adopted therapy.

Multiple sevoflurane exposures don't disturb the T-cell receptor repertoire in infant rhesus monkeys' thymus.

AIMS: Multiple surgical procedures and anesthesia increase the risk of the development in children. However, the influence of such exposures on the developing childhood immunity organs is rarely reported. MATERIALS AND METHODS: High-throughput sequencing of T-cell receptor (TCR) repertoires (TCRseq) from rhesus monkeys' thymus was performed to investigate whether anesthetics could induce de novo antigen recognition via TCR or TCR development impairments. KEY FINDINGS: No significant difference between sevoflurane and control groups regarding VJ gene combinations and diversity of V and J gene was seen, nor was there an obvious change in similar average number of Complementarity Determining Region 3 (CDR3) aa clonotypes. Our analysis of Rank abundance, Gini coefficient, Simpson index, Normalized Shannon Diversity Entropy (NSDE), Morisita-Horn Similarity Index (MHSI) and Bhattacharyya Distance (BD) indicated there is no difference in TCR diversity and similarity. SIGNIFICANCE: These results suggest early events in thymic T cell development and repertoire generation are not abnormality after multiple sevoflurane exposure during childhood. The stabilization of the immune repertoires suggested the safety of sevoflurane in host immune response in children.

Dynamic 3D Locus Organization and Its Drivers Underpin Immunoglobulin Recombination.

A functional adaptive immune system must generate enormously diverse antigen receptor (AgR) repertoires from a limited number of AgR genes, using a common mechanism, V(D)J recombination. The AgR loci are among the largest in the genome, and individual genes must overcome huge spatial and temporal challenges to co-localize with optimum variability. Our understanding of the complex mechanisms involved has increased enormously, due in part to new technologies for high resolution mapping of AgR structure and dynamic movement, underpinning mechanisms, and resulting repertoires. This review will examine these advances using the paradigm of the mouse immunoglobulin heavy chain (Igh) locus. We will discuss the key regulatory elements implicated in Igh locus structure. Recent next generation repertoire sequencing methods have shown that local chromatin state at V genes contribute to recombination efficiency. Next on the multidimensional scale, we will describe imaging studies that provided the first picture of the large-scale dynamic looping and contraction the Igh locus undergoes during recombination. We will discuss chromosome conformation capture (3C)-based technologies that have provided higher resolution pictures of Igh locus structure, including the different models that have evolved. We will consider the key transcription factors (PAX5, YY1, E2A, Ikaros), and architectural factors, CTCF and cohesin, that regulate these processes. Lastly, we will discuss a plethora of recent exciting mechanistic findings. These include Rag recombinase scanning for convergent RSS sequences within DNA loops; identification of Igh loop extrusion, and its putative role in Rag scanning; the roles of CTCF, cohesin and cohesin loading factor, WAPL therein; a new phase separation model for Igh locus compartmentalization. We will draw these together and conclude with some horizon-scanning and unresolved questions.

Identifying the inheritable component of human thymic T cell repertoire generation in monozygous twins.

We have analyzed T cell receptor repertoires in a unique set of thymus samples from a pair of monozygotic twins. While genetics affect the V(D)J rearrangement and generation of junctional sequences, the thymic selections seem largely stochastic and import no detectable inheritable effect at clonal level.

Delineating Human B Cell Precursor Development With Genetically Identified PID Cases as a Model.

B-cell precursors (BCP) arise from hematopoietic stem cells in bone marrow (BM). Identification and characterization of the different BCP subsets has contributed to the understanding of normal B-cell development. BCP first rearrange their immunoglobulin (Ig) heavy chain (IGH) genes to form the pre-B-cell receptor (pre-BCR) complex together with surrogate light chains. Appropriate signaling via this pre-BCR complex is followed by rearrangement of the Ig light chain genes, resulting in the formation, and selection of functional BCR molecules. Consecutive production, expression, and functional selection of the pre-BCR and BCR complexes guide the BCP differentiation process that coincides with corresponding immunophenotypic changes. We studied BCP differentiation in human BM samples from healthy controls and patients with a known genetic defect in V(D)J recombination or pre-BCR signaling to unravel normal immunophenotypic changes and to determine the effect of differentiation blocks caused by the specific genetic defects. Accordingly, we designed a 10-color antibody panel to study human BCP development in BM by flow cytometry, which allows identification of classical preB-I, preB-II, and mature B-cells as defined via BCR-related markers with further characterization by additional markers. We observed heterogeneous phenotypes associated with more than one B-cell maturation pathway, particularly for the preB-I and preB-II stages in which V(D)J recombination takes place, with asynchronous marker expression patterns. Next Generation Sequencing of complete IGH gene rearrangements in sorted BCP subsets unraveled their rearrangement status, indicating that BCP differentiation does not follow a single linear pathway. In conclusion, B-cell development in human BM is not a linear process, but a rather complex network of parallel pathways dictated by V(D)J-recombination-driven checkpoints and pre-BCR/BCR mediated-signaling occurring during B-cell production and selection. It can also be described as asynchronous, because precursor B-cells do not differentiate as full population between the different stages, but rather transit as a continuum, which seems influenced (in part) by V-D-J recombination-driven checkpoints.

Development of Foot-and-Mouth Disease Virus-Neutralizing Monoclonal Antibodies Derived From Plasmablasts of Infected Cattle and Their Germline Gene Usage.

Cattle are susceptible to foot-and-mouth disease virus (FMDV), and neutralizing antibodies are critical for protection against FMDV infection in this species. However, more information is needed on the host specific antigenic structure recognized by the FMDV-specific monoclonal antibodies (mAbs) and on the functional properties of the mAb that are produced in the natural host, cattle. Herein, we characterized 55 plasmablast-derived mAbs from three FMDV-infected cattle and obtained 28 FMDV-neutralizing antibodies by the single B cell antibody technique. The neutralizing mAbs (27/28) mainly recognized conformational epitopes that differ from the well-characterized immunodominant antigenic site 1 of FMDV as defined by murine mAbs. Of these FMDV-neutralizing mAbs, 13 mAbs showed intra-type broadly neutralizing activity against the three topotypes of FMDV serotype O (ME-SA, SEA, and Cathay topotypes). Moreover, all these intra-type broadly neutralizing antibodies competed with sera from FMDV infected or vaccinated cattle, which indicates their binding to native dominant epitopes, as revealed by a blocking ELISA. We further analyzed the germline V(D)J gene usage of the 55 FMDV-specific mAbs and found cattle IgG antibodies containing ultralong HCDR3 were exclusively restricted to usage of the germline gene segment VH 1-7*02. In addition, the restricted germline gene segments of VH 1-7*02 and VL1-47*01 or 1-52*01 pairing were observed in all IgG antibodies with ultralong HCDR3. Furthermore, antibodies with longer HCDR3 were more inclined to display FMDV-neutralizing activity. This study presents a novel method for screening FMDV-specific cattle mAbs which then provide the most useful tools for studying FMDV antigenic structure and variation.

Intra-Vκ Cluster Recombination Shapes the Ig Kappa Locus Repertoire.

During V(D)J recombination, RAG proteins introduce DNA double-strand breaks (DSBs) at recombination signal sequences (RSSs) that contain either 12- or 23-nt spacer regions. Coordinated 12/23 cleavage predicts that DSBs at variable (V) gene segments should equal the level of breakage at joining (J) segments. Contrary to this, here we report abundant RAG-dependent DSBs at multiple Vκ gene segments independent of V-J rearrangement. We find that a large fraction of Vκ gene segments are flanked not only by a bone-fide 12 spacer but also an overlapping, 23-spacer flipped RSS. These compatible pairs of RSSs mediate recombination and deletion inside the Vκ cluster even in the complete absence of Jκ gene segments and support a V(D)J recombination center (RC) independent of the conventional Jκ-centered RC. We propose an improved model of Vκ-Jκ repertoire formation by incorporating these surprisingly frequent, evolutionarily conserved intra-Vκ cluster recombination events.

A Framework for Annotation of Antigen Specificities in High-Throughput T-Cell Repertoire Sequencing Studies.

Recently developed molecular methods allow large-scale profiling of T-cell receptor (TCR) sequences that encode for antigen specificity and immunological memory of these cells. However, it is well-known that the even unperturbed TCR repertoire structure is extremely complex due to the high diversity of TCR rearrangements and multiple biases imprinted by VDJ rearrangement process. The latter gives rise to the phenomenon of "public" TCR clonotypes that can be shared across multiple individuals and non-trivial structure of the TCR similarity network. Here, we outline a framework for TCR sequencing data analysis that can control for these biases in order to infer TCRs that are involved in response to antigens of interest. We apply two previously published methods, ALICE and TCRNET, to detect groups of homologous TCRs that are enriched in samples of interest. Using an example dataset of donors with known HLA haplotype and CMV status, we demonstrate that by applying HLA restriction rules and matching against a database of TCRs with known antigen specificity, it is possible to robustly detect motifs of epitope-specific responses in individual repertoires. We also highlight potential shortcomings of TCR clustering methods and demonstrate that highly expanded TCRs should be individually assessed to get the full picture of antigen-specific response.